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Abstract. The linear optical absorption spectra in neutral boron cluster B6 and cationic Bjt" are calculated 
using a first principles correlated electron approach. The geometries of several low-lying isomers of these 
clusters were optimized at the coupled-cluster singles doubles (CCSD) level of theory. With these optimized 
ground-state geometries, excited states of different isomers were computed using the singles configuration- 
interaction (SCI) approach. The many body wavefunctions of various excited states have been analysed 
and the nature of optical excitation involved are found to be of collective, plasmonic type. 
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1 Introduction 

The cluster science is now a fast emerging field with novel 
properties and tremendous potential for applications. A 
cluster of atoms, or molecules can have few atoms to thou- 
sands of atoms in it, with various structural forms, such 
as spheres, tubes, planar structures etc. Most of the time, 
the properties exhibited by clusters are completely differ- 
ent than their bulk counterpart [1,2]. Boron, in particular, 
is as interesting as carbon, because of its short covalent 
radius and ability to form any structure due to catenation. 
Boron also exhibits structures like nanotubes, fullerenes, 
planar sheets etc. Some of the structures of boron have 
been studied in the context of hydrogen storage [3] . Both 
a and ir aromaticity is observed in many planar boron 
clusters. Some planar boron structures are also found to 
be analogous to hydrocarbons [4]. A circular B^ 9 cluster, 
with a unit of Be wheel in the center behaves as a Wankel 
motor, i.e. the inner B6 wheel rotating opposite to the 
outer B13 ring [5,6]. The all-boron clusters are promising 
candidates as inorganic ligands [7,8]. 

There have been various experimental and theoretical 
investigations of bare boron clusters. Electronic structures 
of boron quasi-planar structures, tubes, wheels and rings 
were studied experimentally using photo-emission spec- 
troscopy by Wang and co-workers [4,9,10,11]. They also 
gave a density functional theory (DFT) based explana- 
tion of their electronic structure. Alexandrova et. al. ex- 
plored the structural and electronic properties and chemi- 
cal bonding in B^ and B6 using anion photoelectron spec- 
troscopy and ah initio calculations [12]. La Placa, Roland 
and Wynne studied abundance spectrum of boron clusters 
generated using laser ablation technique [13]. Ionic boron 
clusters were experimentally studied by Hanley, Whitten 



and Anderson [14]. They studied the fragmentation of big- 
ger clusters into smaller ones by studying the collision- 
induced dissociation of boron clusters. Boustani, on the 
other hand, gave a theoretical description of electronic 
structures of small bare boron clusters using configuration 
interaction method, although with a relatively small Gaus- 
sian basis set [15]. Recently a DFT based study of bare 
boron clusters was done by At is, Ozdogan and Giivenc 
[16]. As far as optical properties of boron clusters are con- 
cerned, there are very few reports available. Marques and 
Botti studied the optical absorption spectra of B 2 o, B 38 , 
B44, Bgo and B92 using time-dependent DFT [17]. To the 
best of our knowledge, no other experimental results are 
available on the optical absorption of boron clusters. 

Our group has recently studied the optical absorption 
in boron clusters B n (n=2 - 5) using a large-scale multi- 
reference configuration interaction method [18]. Since this 
method is quite extensive and scales as N 6 , where N is the 
number of orbitals used in the calculations, it cannot be 
applied for large clusters, or clusters with no symmetry. 
However good insights can be achieved with relatively less 
extensive method known as singles configuration interac- 
tion (SCI), containing only one electron excitations from 
the Hartree-Fock ground state. This method has been ex- 
tensively used for the study of the excited states and opti- 
cal absorption in various other systems [19,20,21,22,23,24]. 
Since optical absorption spectra is very sensitive to the 
structural geometry, the optical absorption spectroscopy 
along with the extensive calculations of optical absorp- 
tion spectra, can be used to distinguish between distinct 
isomers of a cluster. In this report, we present an exten- 
sive calculations of the linear optical absorption spectrum 
of low-lying isomers of B6 and Bg" clusters with different 
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structures. This study along with experimental absorption 
spectra can lead to identification of these distinct isomers. 
Also, in the interpretation of the measured spectra, the 
theoretical insight of the excited states of clusters play an 
important role [25]. 

The remainder of this paper is organised as follows. 
Next section describes the theoretical and computational 
details of the work, followed by section 3, in which re- 
sults are presented and discussed. In the last section we 
present our conclusions, and explore the scope for future 
work. Detailed information regarding the excited states 
contributing to the optical absorption is presented in the 
appendix. 



2 Theoretical and Computational Details 

Different possible arrangements and orientations of atoms 
of the Bq cluster (both neutral and cationic) were ran- 
domly selected for the initial configurational search of ge- 
ometries of isomers. For a given spin multiplicity, the ge- 
ometry optimization was done at a correlated level, i.e. 
at singles doubles coupled cluster (CCSD) level with 6- 
311G(d,p) basis set as implemented in Gaussian 09 [26]. 
Since neutral cluster can have singlet or higher spin mul- 
tiplicity, the optimization was repeated for different spin 
configurations to get the lowest energy isomer. Similarly 
for cationic clusters with odd number of electrons, spin 
multiplicities of 2 and 4 were considered in the optimiza- 
tion. In total, we have obtained 11 neutral Be and 8 cationic 
Bq low- lying isomers. These optimized geometries of neu- 
tral B6 cluster, as shown in Fig. 1, are found to be in good 
agreement with other available reports. Figure 3 shows 
the corresponding geometries of cationic Bq cluster. The 
unique bond lengths, point group symmetry and the elec- 
tronic ground states are given in respective sub-figures. 

For finite systems, such as clusters or quantum dots, 
the ratio J « 1, where a is system size and A is inci- 
dent wavelength of light. In this case, the optical absorp- 
tion cross section of the system corresponding to linear 
absorption can be computed within the electric dipole ap- 
proximation, 

where, a is fine structure constant, Ui n is frequency corre- 
sponding to energy difference of final state \i) and initial 
state \n). f denotes position vector and 7 is linewidth of 
absorption. 

The excited state energies of isomers are obtained us- 
ing the ab initio SCI approach. In this method, differ- 
ent configurations are constructed by replacing one of the 
occupied molecular orbitals in the Hartree-Fock ground 
state by a virtual orbital. Excited states of the system will 
have a linear combination of all such substituted config- 
urations, with corresponding variational coefficients. The 
energies of the excited states will then be obtained by di- 
agonalizing the Hamiltonian in this configurational space 



[27]. The dipole matrix elements are calculated using the 
ground state and the excited state wavefunctions. This is 
subsequently used for calculating the optical absorption 
cross section assuming Lorentzian lineshape, with some 
artificial finite linewidth (cf. Eqn. 1). The contribution of 
wavefunction of the excited states to the absorption peaks 
gives an insight into the nature of optical excitation. 

In one of the earlier reports published elsewhere [18], 
we have extensively studied the dependence of basis sets, 
freezing of Is 2 chemical core and energy cutoff on vir- 
tual orbitals on the computed photoabsorption spectra 
of neutral boron clusters [18]. We have shown that the 
optical absorption spectra of small boron clusters do not 
change even if we freeze the chemical core of boron atoms. 
Also, a careful investigation was done for obtaining the en- 
ergy cutoff of virtual orbitals to be included in the active 
space. Hence, in the present calculations, we have taken 
augmented correlation consistent polarized valence double 
zeta Gaussian basis set [28,29]. We also froze the chemical 
core, and kept only those virtual orbitals with energies of 
1 Hartree or less. 



3 Results and Discussion 

In this section, we discuss the structure and energetics of 
various isomers of neutral and cationic Be cluster, followed 
by discussion of results of computed absorption spectra 
and nature of photo-excitations. 

In the many-particle wavefunction analysis of excited 
states contributing to the various peaks , we have used 
following convention. For doublet systems Hi denotes the 
singly occupied molecular orbital. For triplet systems, two 
singly occupied molecular orbitals are denoted by Hi and 
i^2, while H and L stands for highest occupied molecular 
orbital and lowest unoccupied molecular orbital respec- 
tively. For quartets, the third singly occupied molecular 
orbital is denoted by H3. 



3.1 B 6 

We have found a total of 11 isomers of neutral E>6 cluster 
with stable geometries as shown in the Fig. 1 [30]. The 
relative standings in energy are presented in the Table 1, 
along with point group symmetries and electronic states. 
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Table 1. Point group, electronic state and total energies of different isomers of B6 cluster. 
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Fig. 1. (Color online) Geometry optimized ground state structures of different isomers of neutral B6 clusters, along 
with the point group symmetries obtained at the CCSD level. 
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(a) Planar Ring (Triplet) Isomer 






(b) Incomplete Wheel 
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(c) Bulged wheel Isomer 
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(d) Planar Ring (Singlet) Isomer 
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Fig. 2. (Color online) The linear optical absorption spectrum of Be, calculated using the SCI approach. The many- 
body wavefunctions of excited stated corresponding to the peaks labeled are given in the Appendices A.l-A.ll. For 
plotting the spectrum, a uniform linewidth of 0.1 eV was assumed. 
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The most stable isomer of E>6 cluster has ring-like pla- 
nar structure, with C2/1 point group symmetry. Although 
Be has an even number of electrons, the electronic ground 
state of this isomer is a triplet - an open shell system. 
The equilibrium geometry obtained in our calculation is in 
good agreement with the recently reported values [12,16,31 
The optical absorption spectrum calculated using the SCI 
approach is as shown in the Fig. 2(a). It is mainly charac- 
terized by feeble absortion in the visible range, but much 
stronger absorption at higher energies. The many particle 



wave-functions of excited states contributing to various 
peaks are presented in Table A.l. The first absorption 
peak at 2.84 eV with very low intensity is characterized 
by H — 3 — »• H2 and H — 1 — >• Hi transitions. Due to 
planar nature of the isomer, we can classify the absorp- 
,32^ion into two categories: (a) those with polarization along 
the direction of the plane and (b) polarization perpendic- 
ular to the plane. In this case, it is seen that, in most 
of the cases, the absorption is due to polarizations along 
the plane of the isomer. Also, instead of being dominated 
by single configurations, the wavefunctions of the excited 
states contributing to all the peaks exhibit strong config- 
uration mixing. This is an indicator of plasmonic nature 
of the optical excitations [33] . 

The second low lying isomer of B6 is another planar 
structure resembing an incomplete wheel, i.e. one outer 
atom removed from B7 wheel cluster. This isomer is also 
a triplet system with C2 V symmetry, lying 0.56 eV above 
the global minimum structure. The optimized geometry 
is in good agreement with the other previously available 
reports [31,34]. This is one of the isomers showing opti- 
cal absorption at lower energies (cf. Fig. 2(b)). The many 
particle wave-functions of excited states contributing to 
various peaks are presented in Table A. 2. An intense peak 
at around 2.1 eV is seen in this isomer, characterized by 
Hi —> L + 1 and H 2 — » L transitions. 

A wheel kind of structure, with its center slightly bulged 
out, is found to be the next stable isomer of Bq. A sin- 
glet system with C$ v point group symmetry, lies just 0.83 
eV above in energy as compared to the most stable iso- 
mer. The pentagonal base has bond length of 1.625 A and 
the vertex atom is 1.678 A apart from the corners of the 
pentagon. Other reported values for those bond lengths 
are 1.61 A, 1.66 A [16,34] and 1.61 A, 1.659 A [12,31] re- 
spectively. The optical absorption spectrum is presented 
in Fig. 2(c). The many-particle wave functions of excited 
states contributing to various peaks are presented in Table 
A. 3. The onset of the spectrum occurs near 3.76 eV, with 
polarization in the plane of the pentagonal base, charac- 
terized by excitations H — 1 — » L + 6 and H — » L + 6 with 
equal contribution. 

A planar ring like structure, resembling the global min- 
imum one; however, with singlet state and C s symmetry, 
is the next low lying isomer of Bq cluster. The optimized 
geometry is in good agreement with Ref. [34] . The absorp- 
tion spectrum is presented in Fig. 2(d) and corresponding 
many-particle wavefunctions of various excited states are 
presented in Table A. 4. The absorption starts at the end of 
visible range, characterized by H—l —> L and H — >• L+23. 

An octahedron structure with Oh point group symme- 
try is the next stable isomer of neutral B6. Each side of the 
octahedron is found to be 1.7 A as compared to the 1.68 
A reported by [16] and [34]. The many-particle wave func- 
tions of the excited states correponding to various peaks 
(cf. Fig. 2(e)) are presented in Table A. 5. A very feeble 
absorption at 3.6 eV opens the spectrum, mainly charac- 
terized by excitations H — 1 — >> L and H — 1 — >> L + 1 with 
equal contribution. 
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Next isomer is previously unreported, with structure 
of a saddle threaded with dimer from top. It lies just 0.04 
eV above the previous octahedron isomer. However the 
optical absorption spectrum (cf. Fig. 2(f)) is completely 
different. A narrow energy range hosts all the peaks. The 
onset of spectrum occurs near 2.7 eV, with major con- 
tribution from H — 1 H 2 and Hi — >• L + 9 (cf. Table 
A.6). 

An isosceles triangle connected to a linear chain of 
boron atoms forms the next isomer. This structure with 
C2V symmetry and a triplet electronic state have been 
reported in [34], which is in close agreement with our re- 
sults. The optical absorption spectrum (cf. Fig. 2(g)) has 
distinctive closly lying peaks at 5.42 eV and 5.66 eV. The 
many-particle wavefunctions of excited states correspond- 
ing to various peaks are presented in Table A. 7. Config- 
urations H — 3 — >• L and H — 1 — >• L + 1 contribute pre- 
dominently to those closly lying peaks respectively. 

A structure with two out of plane isosceles triangles 
joined together is found to be one of the isomers. The ge- 
ometry has isosceles triangle with lengths 1.62 A, 1.62 A 
and 1.53 A, while two such triangles are joined by a bond 
of length 1.66 A. The respective numbers reported by Ref. 
[34] are 1.60 A, 1.60 A, 1.50 A and 1.647 A respectively. 
The optical absorption spectrum contains many low inten- 
sity peaks except for strongest one at 5.87 eV, as presented 
in Fig. 2(h). The onset of the spectrum occurs at 1 eV - 
a peak equally dominated by H — >> L and H — 1 — >• L (cf. 
Table A.8). 

An almost degenerate structure forms the next isomer, 
lying just 0.009 eV above the previous isomer. Contrary 
to the previous one, this geometry is completely planar 
and is a triplet system, with C2 V point group symmetry. 
Probably because of such a strong near-degeneracy, this 
isomer has not been reported in the literature before. The 
many-particle wave functions of the excited states corre- 
ponding to various peaks (cf. Fig. 2(i)) are presented in 
Table A. 9. The spectrum opens with a very feeble peak 
H — >• L and H — l^L + las dominant contributions to 
the excited state wavefunction. 

Convex bowl shaped isomer and a perfect linear chain 
are found very high in energy, ruling out their existence 
at room temperature. The optical spectra are presented in 
Figs. 2(j) and 2(k) respectively. The corresponding many- 
particle wavefunctions of excited states of various peaks 
are presented in Table A. 10 and A. 11. The onset of ab- 
sorption spectrum of convex bowl isomer occurs at 1.58 
eV with H — 1 — >• L + 1 as the most dominant configura- 
tion. The bulk of the oscillator strength of the spectrum of 
linear isomer is carried by H — 1 — >• L + 3 and H — >> L + 2 
having equal contributions. 



3.2 B+ 

We have found a total of 8 isomers of cationic (Bg") clus- 
ter with stable geometries as shown in the Fig. 3. The 
relative standings in energy are presented in the Table 2, 
along with point group symmetries and electronic states. 



Table 2. Point group, electronic state and total energies 
of different isomers of Bg" cluster. 
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Cationic clusters show activity in the visible range, con- 
trary to their neutral counterpart. Also, in most of the 
cases the geometry of the neutral isomer is retained, re- 
flected in the fact that some peaks show up in the optical 
absorption spectra at the same energies as those in the 
neutral cluster. 

The most stable isomer of Bg" cluster is a planar ring- 
type of structure, with C s point group symmetry. This is 
in contrast to the other reported geometries which have 
symmetry [31,32,34]. A slight difference in the orien- 
tation makes it less symmetric. However, the bonds lengths 
obtained are in good agreement with those with T)2h sym- 
metric geometry cited above. The optical absorption spec- 
trum calculated using SCI approach is as shown in the 
Fig. 4(a) and corresponding many particle wavefunctions 
of excited states contributing to the various peaks are pre- 
sented in Table A. 12. Similar to the neutral counterpart, 
this isomer also has very feeble absorption in the visible 
range, with polarization perpendicular to the plane of the 
isomer. The dominant contribution to the peak at 0.63 eV 
comes from Hi —> L and H — 2 ^ L + 16 . The spectrum 
is slightly red-shifted with respect to the neutral counter- 
part. 

Bulged wheel structure is the next low lying isomer of 
Bg" with just 0.023 eV above the global minimum. How- 
ever, as compared to the neutral one, this geometry has Ci 
symmetry due to the significant bond length reordering. 
Our results are consistent with the results of Refs. [31,34]. 
The optical absorption spectrum is presented in Fig. 4(b). 
The many-particle wave functions of excited states con- 
tributing to various peaks are presented in Table A. 13. 
The spectrum is distinctly different with bulk of the os- 
cillator strength carried by peaks at 6 eV. The onset of 
spectrum occurs at 1.67 eV dominated by H — 2 — >> Hi 
and Hi — >> L + 3 configurations. 

Another ring-like structure with T)2h symmetry and 
doublet multiplicity lies next in the energy order. The 
six faced benzene type structure has 1.66 A and 1.53 A 
as unique bond lengths, which are somewhat larger than 
those reported in the literature [31, 34]. The optical absorp- 
tion spectrum is presented in Fig. 4(c). The many-particle 
wave functions of excited states contributing to various 
peaks are presented in Table A. 14. All absorption peaks 
are due to the polarization along the plane of the isomer. 
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(a)C„ 2 A" (b)d, 2 A (c)D 2h , 2 A g (d) C 2v , 4 A 2 

Planar ring (I) Bulged wheel Planar ring (II) Incomplete wheel 




(h) D 2 h, 2 B2 9 Planar trimers 



Fig. 3. (Color online) Geometry optimized ground state structures of different isomers of cationic Bg" clusters, along 
with the point group symmetries obtained at the CCSD level. 



The onset of spectrum occurs at 1.2 eV with a very feeble 
peak with dominant contributions from H — 1 —> L and 
H-2^H 1 . 

Next low lying isomer is a planar incomplete wheel 
structure with Q 2 v point group symmetry and quartet 
multiplicity. This is consistent with results of Ref. [34]. 
The optical absorption (cf. Fig. 4(d)) starts at 2.94 eV, 
with polarization transverse to the plane of the isomer. 
The many particle wave-functions of excited states con- 
tributing to various peaks are presented in Table A. 15. 
The configurations contributing to the first peak are Hi — >> 
L + 1 and H 2 — » L, where Hi and H 2 are singly occupied 
molecular orbitals of the quartet system. 

Another isomer with the same point group symmetry 
and multiplicity as that of the previous one, but having 
a geometry of linear chain with an isosceles triangle at 
the end, is the next low lying isomer of cationic Bg". Our 
results about geometry are in good agreement with the 
Ref. [34]. The optical absorption spectrum (cf. Fig. 4(e)) 
has very few peaks. The onset of spectrum occurs at 3.2 
eV, with polarization along the plane of the isomer, and 
has dominant contribution from H — 3 —> H3 and H — 3 —> 
L (cf Table A. 16). 

Tetragonal bipyramid forms the next stable isomer of 
cationic Bg", with point group symmetry and doublet 
multiplicity. This is in good agreement with the geome- 
tries reported in Refs. [31] and [34]. The optical absorp- 
tion spectrum (cf Fig. 4(f)) has peaks in the visible range 
at 1.58 eV and 1.97 eV, with dominant contribution from 



Hi L + 2 and H - 1 -> L + 2 respectively (cf. Table 
A.17). 

Two more structures were found stable i.e. (a) a pla- 
nar structure with two trimers joined together and, (b) 
a linear one. These isomers are much above the global 
minimum energy, it rules out their room temperature ex- 
istence. In the linear isomer the absorption begins with a 
small peak at 0.93 eV, with dominant contribution from 
H — 3 — » H3 and H — 1 — »> L configurations. In case of 
planar trimers structure, the spectrum seems to be red 
shifted as compared to the neutral one. The first peak is 
found at 0.77 eV with H - 2 Hi and Hi L + 4 as 
dominant contribution to the wavefunction of the excited 
state. 

4 Conclusion and Outlook 

A large number of randomly selected initial structures of 
neutral Be and cationic Bg" clusters are taken into consid- 
eration for locating the global and local minimas on the 
potential energy curves. A careful geometry optimization 
is done for all those structures at a correlated level. The 
optical absorption spectra of different low-lying isomers 
of both neutral and cationic isomers are reported here. 
A singles configuration interaction approach was used to 
compute excited state energies and the absorption spec- 
tra of various clusters. Spectra of cationic clusters appear 
slightly red-shifted with respect to the neutral one. Differ- 
ent isomers exhibit distinct optical response, even though 



8 Ravindra Shinde and Alok Shukla: Optical absorption in boron clusters B6 and B^~: A first principles CI approach 




02468 10 02468 10 02468 10 



(a) Planar Ring (I) Isomer (b) Bulged wheel Isomer (c) Planar Ring (II) Isomer 




2 4 6 S 10 02468 10 02468 10 



(d) Incomplete Wheel (e) Threaded trimer (f) Tetragonal bipyramid 




Energy (eV) Energy (eV) 
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Fig. 4. (Color online) The linear optical absorption spectrum of Bg", calculated using the SCI approach. The many- 
body wavefunctions of excited stated corresponding to the peaks labeled are given in the Appendices A. 12- A. 19. For 
plotting the spectrum, a uniform linewidth of 0.1 eV was assumed. 



they are isoelectronic and many of them are almost degen- 
erate. This signals a strong-structure property relation- 
ship, which can be exploited for experimental identifica- 
tion of these isomers; something which is not possible with 
the conventional mass spectrometry. A strong mixture of 
configurations in the many-body wave functions of vari- 
ous excited states are observed, indicating the plasmonic 
nature of the photoexcited states [33]. 



Since aluminum also has the same number of valence 
electrons, it will be interesting to compute its optical ab- 
sorption spectra, and compare them with those of boron 
cluster with the help of many-body wavefunctional anal- 
ysis. The results of such calculations done by us will be 
communicated in near future [35]. 



The author would like to acknowledge the Council of Scien- 
tific and Industrial Research (CSIR), India, for their financial 
support (SRF award No. 09/087/(0600)2010-EMR-I and travel 
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A Excited State CI Wavefunctions, Energies 
and Oscillator Strengths 

In the following tables, we have given the excitation energies 
(with respect to the ground state), and the many body wave- 
functions of the excited states, corresponding to the peaks in 
the photoabsorption spectra of various isomers listed in Fig. 
1 and Fig. 3 , along with the oscillator strength /i 2 of the 
transitions, 

fi2 = l^(E 2 - E 1 )Y,\(m\d i \G)\ 2 (2) 
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where, \m) denotes the excited state in question, \G), the ground 
state, and di is the i-th Cartesian component of the electric 
dipole operator. The single excitations are with respect to the 
reference state as given in respective tables. 
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Table A.l. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of - planar ring (triplet) isomer (cf. Fig. 2(a)). The subscript || in the peak number 
denotes the absorption due to light polarized in the plane of isomer. In the wave function, the bracketed numbers are 
the CI coefficients of a given electronic configuration. Symbols #i,i?2 denote SOMOs discussed earlier, and H, and 
L, denote HOMO and LUMO orbit als respectively. Note that, the reference state does not correspond to any peak, 
instead it represents the reference state from which singles excitations are occuring. 
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H - 


1 - 


^ L + 8)(0.5137) 


IV|| 


6.62 


0.9942 


H - 


1 - 


-> L + 8) (0.5440) 








H - 


1 - 


-> L + 3) (0.4900) 


V|| 


7.96 


0.5799 


H- 


1 - 


^ L + 3) (0.6800) 








H - 


1 - 


^ L + 8) (0.5298) 



Table A. 2. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of B6 - incomplete wheel isomer (cf. Fig. 2(b)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of wheel base, respectively. 
The rest of the information is the same as given in the caption for Table A.l. 



Peak 


E (eV) 


fi2 Wave Function 


Reference 






HhHl) 


III 


2.09 


0.1460 


Hi L + l) (0.8102) 
H 2 -+ L)(0.2871) 


H|| 


3.28 


0.0406 


H 2 L)(0.9011) 
H- 1 -> Hi) (0.2359) 


III± 


3.86 


0.0247 


H - 3 -+H 2 ) (0.9301) 
if -3^ L + 27) (0.2154) 


IV|| 


4.28 


0.0849 


H - 2 ->H 2 ) (0.7926) 
if - 1 Hi) (0.3875) 


V|| 


5.26 


0.1008 


H- 5 ->H!) (0.4241) 
if - 1 Hi) (0.4168) 


Vly 


5.56 


0.2772 


Hi + 16) (0.5724) 
Hi -+ L + 5) (0.3863) 


VII|| 


5.96 


0.0677 


Hi + 18) (0.5153) 
Hi -> L + 8) (0.3795) 


VIII|| 


6.57 


0.1286 


H-l-> L) (0.5208) 
H - 5 ->L) (0.5092) 
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Table A. 3. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Be - bulged wheel isomer (cf. Fig. 2(c)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of wheel base, respectively. In 
the wave function, the bracketed numbers are the CI coefficients of a given electronic configuration. Symbols H and 
L denote HOMO and LUMO orbitals respectively. Excitations are with respect to Hartree Fock reference state. 



Peak E (eV) f\2 Wave Function 



1,1 3.76 


0.2625 


\H- 


1 - 


+ L + 6) (0.5622) 






\H 


> L + 6) (0.5622) 






11 


1 - 


+ L + 2) (0.4459) 






\H - 


> L + 2) (0.4459) 


II II 4.38 


0.0657 


\H- 


1 - 


^ L + 3) (0.4998) 






H — 
1 1 


> L + 3) (0.4998) 






\H- 


1 - 


+ L + 4) (0.4981) 






\H - 


> L + 4) (0.4981) 


IIii, 5.07 


2.1039 


\H- 


2 - 


-> L + 1) (0.7653) 






\H- 


2 - 


+ L) (0.7653) 






\H- 


2 - 


+ L + 8) (0.2857) 






\H - 


2 - 


-> L + 7) (0.2857) 


IV± 5.71 


0.0810 


\H- 


2 - 


-> L + 6) (0.6266) 






\H- 


2 - 


+ L + 2) (0.5564) 


Vy 7.05 


1.0550 


\H - 


> L + 5) (0.7645) 






\H- 


1 - 


+ L + 5) (0.7645) 






\H - 


> L + 6) (0.4763) 






\H- 


1 - 


+ L + 6) (0.4763) 



Table A. 4. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of B6 - planar ring (singlet) isomer [cf. Fig. 2(d)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, respectively. 
The rest of the information is the same as given in the caption for Table A. 3. 



Peak E (eV) f\2 Wave Function 



III 


3.29 


0.5397 


H - 
H - 


1 -> L) (0.9174) 
> L + 23) (0.2504) 


Illl 


5.46 


1.0153 


H - 
H - 


> L + 7) (0.4991) 

> L + 23) (0.4093) 


111,1 


5.77 


0.1380 


H — 
H - 


1 L + 2) (0.5479) 
1 L + 8)(0.5181) 




6.04 


0.8199 


H - 
H - 


> L + 7) (0.7472) 

> L + 12)(0.3140) 


V|| 


6.93 


2.2212 


H - 
H - 


> L + 16) (0.4623) 
y L + 12) (0.4600) 


VI|| 


7.31 


0.6003 


H - 
H - 


> L + 16) (0.6363) 
y L + 12) (0.5687) 
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Table A. 5. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of Be - octahedron isomer (cf. Fig. 2(e)). The subscripts || and _L, in the peak number 
denote the absorption due to light polarized in, and perpendicular to the plane of pyramidal base, respectively. The 
rest of the information is the same as given in the caption for Table A.l. 



Peak E (eV) f\2 Wave Function 



Reference 






\HhHl) 






3.63 


0.0047 


\H - 1 -> 
if -1 -> 
if -2 -> 
if -2-> 


L)(0.8716) 
L + l)(0.8716) 
L + 1) (0.2638) 
L) (0.2638) 




5.56 


0.0129 


|if - 1 -> 
if -1 -> 
if- 1 -> 
if- 1 -> 


L + 1) (0.5440) 
L) (0.5440) 
L) (0.5277) 
L + 1) (0.5277) 




6.05 


0.0357 


|if -1 -> 
if -1 -> 
if- 1 -> 
if -1 -> 


L) (0.6603) 
L + 1) (0.6603) 
L + 1) (0.4383) 
L) (0.4383) 


IV||,± 


7.65 


0.0149 


|if -1 -> 
if- 1 -> 
if -2 -> 
jif -2 -> 


L)(0.6741) 
L + 1) (0.4383) 
if 2 ) (0.5346) 
L + 1) (0.5346) 



Table A. 6. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of B 6 - threaded tetramer isomer (cf. Fig. 2(f)). The subscript x in the peak number 
denote the absorption due to light polarized along the long axis, and, y, z denotes polarization perpendicular to it. 
The rest of the information is the same as given in the caption for Table A.l. 



Peak 


E (eV) 


fl2 


Wave Function 


Reference 






\HbHl) 


I* 


2.74 


0.1116 


|if - 1 if 2 ) (0.5535) 
ifi -> L + 9) (0.3854) 


Hx 


3.33 


0.2917 


|if - 1 if 2 ) (0.7564) 
ifi L + 9) (0.2492) 
ifi L + 3) (0.2465) 


III* 


3.90 


0.0480 


|ifi L) (0.8447) 

if - 1 L + 9)(0.2022) 


IV, 


4.33 


0.1507 


|if - 1 L + 9)(0.4175) 
if - 1 -> L + 3) (0.4043) 


V, 


5.09 


0.0338 


| if -3-> L + 9)(0.5013) 
if-3-> L + 3) (0.4697) 


VI, 


5.43 


0.0852 


|if - 1 -> L + 1)(0.7137) 
if -2 ^L) (0.2694) 
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Table A. 7. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of Be - threaded trimer isomer (cf. Fig. 2(g)). The subscript || , in the peak number 
denotes the absorption due to light polarized along the long axis of the isomer. The rest of the information is the same 
as given in the caption for Table A.l. 



Peak E (eV) f\2 Wave Function 



Reference 






\HbHl) 




III 


5.42 


1.4376 


\H-3^ 
\H- 1 -> 


L) (0.8262) 
L + 1)(0.3212) 


H|| 


5.66 


2.6846 


\H- 1 -> 
tf-3-> 


L + 1)(0.6710) 
L) (0.3680) 


III,, 


6.25 


0.5077 


\H -4^ 
\H- 1 -» 
\H- 1 -> 


L + 1) (0.6263) 
L + 11) (0.2702) 
L + l) (0.2701) 
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Table A. 8. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of - twisted trimers isomer (cf. Fig. 2(h)). The subscripts || and _L, in the 
peak number denote the absorption due to light polarized along, and perpendicular to the long axis of the isomer, 
respectively. The rest of the information is the same as given in the caption for Table A. 3. 

Peak E (eV) fi2 Wave Function 



Iy 1.02 0.0422 

II ± 2.22 0.1279 

Hill 3.58 0.1870 

IVy 4.72 0.4337 

V± 5.26 0.3570 

VI_l 5.87 3.4368 

VII|| 6.39 0.2265 

VIII± 6.98 0.6645 

IX|| 7.25 0.3140 



H - 


> L) (0.6550) 


H - 


1 - 


* L) (0.6550) 


H - 


>L + 1)(0.5612) 


H - 


1 - 


-> L + 1)(0.5612) 


H - 


2 - 


L) (0.7650) 


H - 


3 - 


* L + l)(0.5316) 


H - 


2 - 


^ L + 2) (0.4374) 


H - 


5 - 


L) (0.4374) 


H - 


4 - 


* L) (0.4371) 


H - 


2 - 


^ L + 3) (0.4371) 


H - 


2 - 


^ L + 2) (0.4874) 


H - 


2 - 


L + 3) (0.4874) 


H - 


> L + 26) (0.3268) 


H - 


1 - 


^ L + 26) (0.3268) 


H - 


5 - 


* L + 3)(0.5117) 


H - 


4 - 


L + 2)(0.5010) 


H - 


> L + 3) (0.5962) 


H - 


1 - 


-> L + 2) (0.5962) 


H - 


3 - 


L + 2) (0.4908) 


H - 


3 - 


-> L + 3) (0.4908) 


H - 


2 - 


L + 3) (0.4561) 


H - 


2 - 


^ L + 2) (0.4561) 


H - 


3 - 


L + 1) (0.6479) 


H - 


2 - 


-> L) (0.5508) 


H - 


> L + 4) (0.5209) 


H - 


1 - 


L + 4) (0.5209) 


H - 


> L + 7) (0.4072) 


H - 


1 - 


^ L + 7) (0.4072) 



Ravindra Shinde and Alok Shukla: Optical absorption in boron clusters E>6 and B^~: A first principles CI approach 15 

Table A. 9. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg - planar trimers isomer (cf. Fig. 2(i)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, respectively. 
The rest of the information is the same as given in the caption for Table A. 3. 

Peak E (eV) fi2 Wave Function 



I_L 


0.97 


0.0393 




-> L) (0.7804) 








\H 


- 1 - 


^ L + 1) (0.5529) 


II 


2.22 


0.1261 


\H 


-2 - 


* L) (0.7603) 








\H 


-4 - 


+ L + 1) (0.5634) 


Illy 


3.57 


0.1293 


\H 


-2 - 


-> L + 2) (0.6070) 








\H 


-3 - 


+ L) (0.5440) 


IV|| 


4.67 


0.7257 


\H 


-3 - 


+ L + 2) (0.7242) 








\H 


L + 12) (0.3824) 








\H 


-2 - 


^ L + 2) (0.5131) 








\H 


-3 - 


+ L) (0.4653) 


V|| 


6.43 


4.5721 


\H 


-> L + 12) (0.5256) 








\H 


-3 - 


^ L + 2) (0.4859) 


VI|, 


6.99 


1.8466 


\H 


-4 - 


+ L + 1) (0.5794) 








\H 


L + 7)(0.5156) 


VII|| 


7.34 


0.6386 


\H 


L + 12) (0.6754) 








\H 


-+ L + 7) (0.4922) 
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Table A. 10. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks in 
the linear absorption spectrum of - convex bowl isomer (cf. Fig. 2(j)). The subscripts || and _L, in the peak number 
denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, respectively. The rest 
of the information is the same as given in the caption for Table A. 3. 



Peak 


E (eV) 


fi2 Wave Function 


III 


1.58 


0.0486 


H- 


1 -> L + 1) (0.9774) 




1.80 


0.0679 


H - 
H - 


> L + 2) (0.9079) 
1 L + 3) (0.2748) 


Ill,, 

II 


2.43 


0.6023 


H - 
H - 


> L) (0.8644) 
4 L + 1) (0.4364) 


IVii 

II 


2.89 


0.1811 


H - 
H - 


1 L) (0.8058) 
3 L + 1) (0.5421) 


Vn 
v II 


4.09 


0.1602 


H - 
H - 


4 L) (0.5608) 
3 -> L + 1)(0.5347) 


VI||,± 


5.13 


0.0649 


H - 
H - 


4 + 3) (0.9070) 
3 L + 2)(0.1753) 




6.21 


1.4363 


H - 
H - 


> L + 7>(0.4118) 

> L + 4) (0.3622) 


VIII_L 


6.39 


2.2019 


H - 
H - 


> L + 6) (0.4335) 
2 -> L + 3)(0.4081) 


IX,, 


6.90 


0.1476 


H - 
H - 


> L + 8) (0.4435) 

> L + 13) (0.3534) 



Table A. 11. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of B6 - linear isomer {cf. Fig. 2(k)). The subscripts || and _L, in the peak number 
denote the absorption due to light polarized along, and perpendicular to the axis of the isomer, respectively. The rest 
of the information is the same as given in the caption for Table A. 3. 



Peak E (eV) f\2 Wave Function 



III 


5.51 12.8358 \H - 


1 L + 3)(0.6510) 




\H~ 


> L + 2) (0.6489) 


II± 


6.51 1.6532 \H - 


> L + 4) (0.7148) 




\H~ 


1 -+ L + 4) (0.7148) 




\H- 


3 L + 8) (0.3847) 




\H~ 


2 L + 8) (0.3847) 
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Table A. 12. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg" - planar ring isomer (cf. Fig. 4(a)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, respectively. In 
the wave function, the bracketed numbers are the CI coefficients of a given electronic configuration. Symbols Hi 
and L denote HOMO, SOMO and LUMO orbitals respectively. 



Peak 


E (eV) 


/12 Wave Function 


Reference 










l± 


0.63 


0.0075 


Hx - 
H - 


-> L) (0.9667) 
2 -> L + 16)(0.1229) 


nil 


4.43 


0.9637 


H- 
H! - 


1 L) (0.7448) 
^ L + 10) (0.4227) 


in,. 


5.18 


0.0565 


H - 
H- 


2 -+ L + 1) (0.6283) 
1 -> L + 1) (0.5264) 


IV|, 


5.64 


0.4310 


H - 
H- 


1 L + 2) (0.6350) 
3 L + 3) (0.5173) 


V|| 


6.06 


0.1576 


H- 
H - 


3 L + 3) (0.5624) 
1 -+ L + 3) (0.4278) 
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Table A. 13. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - bulged wheel isomer (c/. Fig. 4(b)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of the wheel, respectively. The 
rest of the information is the same as given in the caption for Table A. 12. 



Peak 


E (eV) 


fl2 


Wave Function 


Reference 






\h\) 




III 

ii 


1.67 


0.0159 


\H- 
\Hi - 


2 Hi) (0.9568) 
-> L + 3) (0.1270) 


n ll 

ii 


3.25 


0.0061 


\m - 

\H! - 


-> L) (0.5877) 
L + 3) (0.3945) 


III,, 

ii 


3.65 


0.0021 


\H! - 

\Hi ~ 


-> L + 2)(0.8106) 
-> L + 1) (0.2830) 


IVy 


4.20 


0.0998 


\H- 
\H - 


1 L) (0.5760) 

2 L) (0.5743) 


V|| 


4.46 


0.0312 


\H- 
\H X - 


1 -> L + 2) (0.7945) 
-> L)(0.2711) 


VI|, 


4.69 


0.0417 


\H- 
\H - 


1 L + 2) (0.7945) 

2 -> L)(0.2711) 


VII|| 


4.91 


0.0453 


\H- 
\H 


2 L + 1)(0.6517) 
1 -> L + 1) (0.4364) 


VIII|, 


5.60 


0.4219 


11 

\H- 


1 L + 3) (0.6402) 

2 -> L) (0.4002) 


IX||,± 


5.77 


0.2700 


\H- 
\H- 
\H- 
\H - 


2 L + 2) (0.8807) 

1 L + 4)(0.6632) 

2 L + 1) (0.3734) 
2 -> L) (0.2952) 


X|| 


6.00 


0.5544 


\H- 
\H- 


3 L) (0.4788) 
2 L + 1)(0.4263) 


XI,, 


6.27 


0.3573 


\H- 
\H- 


2 -> L + 3) (0.6862) 
4 L) (0.4386) 
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Table A. 14. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - planar ring (II) isomer (cf Fig. 4(c)). The subscript || in the peak number 
denote the absorption due to light polarized in the plane of the isomer. The rest of the information is the same as 
given in the caption for Table A. 12. 



Peak E (eV) f\2 Wave Function 



Reference 






tfi 1 ) 






III 


1.18 


0.0182 


H - 


1 


L) (0.6957) 








H - 


2 


iZi) (0.6943) 


H|| 


3.47 


1.5345 


H - 


1 


-> L) (0.7424) 








H - 


2 


ifi) (0.4272) 


III,, 


4.70 


0.1042 


if - 


2 


L + 1) (0.7037) 








if - 


4 


L) (0.4979) 


IV|| 


4.93 


0.1114 


if - 


1 


L + 3) (0.9077) 








if - 


5 


L) (0.1606) 


V|| 


5.12 


0.0673 


if - 


4 


-> L) (0.7687) 








if - 


3 


-> L + 2)(0.5208) 


VI|| 


6.43 


0.8286 


ifi - 




L + 6) (0.6655) 








if - 


1 


L + 3) (0.4625) 



Table A. 15. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - incomplete wheel (quartet) isomer (cf Fig. 4(d)). The subscripts || and _L, 
in the peak number denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, 
respectively. In the wave function, the bracketed numbers are the CI coefficients of a given electronic configuration. 
Symbols if, H u if 2 , H 3 and L denote HOMO, SOMOs and LUMO orbitals respectively. 

Peak E (eV) fi2 Wave Function 
Reference \Hl; if^; H\) 



I± 


2.94 


0.0608 


ifi - 
H2 - 


-+ L + l)(0.7513) 
L) (0.6261) 


nil 


3.25 


0.1207 


if - 
if - 


5 if 3 ) (0.7994) 
4 — >• L) (0.3377) 


HI_l 


4.84 


0.0248 


if 3 - 
ifi - 


-> L) (0.8592) 
-> L + 1) (0.2632) 


IV|| 


5.60 


0.0552 


if - 
ifi - 


3 ->> #2) (0.6397) 
-+ L + 3) (0.4131) 


V|| 


6.8 


0.0674 


if - 
if - 


3 ifi) (0.6876) 

4 -> L) (0.4045) 


VI|, 


7.25 


0.1363 


if - 
if 2 - 


5 if 2 ) (0.7232) 
-> L) (0.2883) 


VII± 


7.39 


0.0336 


if - 
ifi - 


3 L)(0.4522) 
-+ L + 5) (0.3485) 
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Table A. 16. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - threaded trimer (quartet) isomer (cf. Fig. 4(e)). The subscripts || and _L, in 
the peak number denote the absorption due to light polarized along, and perpendicular to the long axis of the isomer, 
respectively. The rest of the information is the same as given in the caption for Table A. 15. 



Peak 


E (eV) 


/l2 


Wave Function 


Reference 






\Hi]HI]HI) 


III 


3.18 


0.4369 


\H - 3 -+H 3 ) (0.6993) 
\H- 3 (0.5802) 




4.48 


0.1117 


\H 3 -> L + 2) (0.8855) 
\H- 3 ->L) (0.7012) 
if-4^L + 2)(0.4816) 
#-3^L + 3)(0.3931) 


III,, 


5.03 


2.3151 


\H 3 -> L) (0.6337) 
\H 2 -> L + 1) (0.5640) 



Table A. 17. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - tetragonal bipyramid isomer (cf. Fig. 4(f)). The subscripts || and _L, in the 
peak number denote the absorption due to light polarized in, and perpendicular to the square plane of bipyramid, 
respectively. The rest of the information is the same as given in the caption for Table A. 12. 



Peak 


E (eV) 


fl2 


Wave Function 


Reference 






l^i 1 ) 






l± 


1.58 


0.2525 


|#i - 


->> 


L + 2) (0.7189) 








\H- 


1 


-> L + 2)(0.6154) 


H|| 


1.97 


0.0803 


\H! - 




L + 1) (0.7492) 








\H- 


1 


-> L + 1)(0.5593) 




3.87 


0.0272 


\H- 


2 


-> L) (0.8717) 








\H- 


1 


L) (0.8684) 








\H- 


2 


L + 1) (0.4082) 








\H! - 




L) (0.3497) 


IV|, 


4.36 


0.0287 


\H- 


2 


-> L + 2)(0.8419) 








\H- 


1 


L) (0.3385) 


V|| 


4.55 


0.0440 


\H- 


2 


L + 2)(0.9001) 








\H- 


4 


-> L + 2) (0.8807) 








\H- 


1 


L) (0.2555) 








\H- 


5 


L) (0.1997) 


VI II,^ 


5.04 


0.8028 


\H- 


4 


L + 2)(0.8991) 








\H - 


1 


L + 2)(0.6152) 








\H- 


4 


L + 3) (0.5964) 








\H- 


1 


L) (0.5334) 








\H- 


2 


L + 2) (0.4299) 








\H- 


1 


L + 3) (0.3269) 
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Table A. 18. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg" - linear (quartet) isomer (c/. Fig. 4(g)). The subscripts || and _L, in the 
peak number denote the absorption due to light polarized along, and perpendicular to the long axis of the isomer, 
respectively. The rest of the information is the same as given in the caption for Table A. 15. 



Peak 


E (eV) 


/12 


Wave Function 


Reference 






\H\ ; Hi ; Hi) 


III 


0.93 


0.1533 


\H - 3 -+H 3 ) (0.8624) 
\H 2 -> L)(0.3187) 


II± 


3.62 


0.0677 


\H 2 L) (0.7478) 
\H! -> L) (0.6474) 


III,, 


5.02 


0.6861 


\H- 3 + (0.9263) 
if- 3 + 12) (0.2264) 


IV|| 


5.60 


9.8432 


\H-4->L) (0.6486) 
\H- 3 + (0.3521) 



Table A. 19. Excitation energies, E, and many-particle wave functions of excited states corresponding to the peaks 
in the linear absorption spectrum of Bg~ - planar trimers isomer (cf. Fig. 4(h)). The subscripts || and _L, in the peak 
number denote the absorption due to light polarized in, and perpendicular to the plane of the isomer, respectively. 
The rest of the information is the same as given in the caption for Table A. 12. 



Peak 


E (eV) 


fi2 Wave Function 


Reference 






Hi) 




III 


0.77 


0.2915 


H - 
#i - 


2 -> Hi) (0.9647) 
L + 4)(0.1605) 


II± 


1.33 


0.0108 


#i - 
H - 


-> L) (0.8340) 
2 L + 1) (0.5093) 


IH± 


2.08 


0.0378 


H- 
H! - 


2 L + 1) (0.9744) 
-+ L)(0.1888) 


IV_L 


2.51 


0.0103 


H - 
H - 


3 -> Hi) (0.9414) 
1 -> L + 4)(0.2905) 


V|| 


2.86 


0.1131 


H - 
H - 


1 L) (0.7701) 
3 L + 1)(0.5892) 


VI|| 


4.16 


0.1267 


H- 
H - 


1 L + 2)(0.6627) 
4 L) (0.5105) 


Vlly 


5.30 


0.4383 


Hi - 
H - 


^ L + 9) (0.5023) 
1 L + 2) (0.4505) 
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